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Abstract 

Residual  thermal  stresses  in  the  electrolyte  of  the  anode-supported  planar  solid  oxide  fuel  cells  (SOFCs)  were  estimated  experimentally 
and  numerically.  The  X-ray  diffraction  method  was  employed  to  measure  the  residual  stress  in  the  cell.  To  estimate  the  residual  stress 
precisely,  a  synchrotron  radiation  was  used  as  the  best  means  of  X-ray  source.  The  wavelength  of  the  radiation  beam  used  for  the  stress 
measurements  was  0.154nm,  and  the  diffraction  peak  of  YSZ  (5  3  1)  plane  was  used  for  the  stress  measurement.  The  estimated  residual 
stresses  in  the  electrolyte  were  compressive  stresses  of  around  650  MPa  at  room  temperature.  Numerical  calculation  for  the  residual  stress 
was  also  carried  out. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Anode-supported  solid  oxide  fuel  cells  (SOFCs)  are  suit¬ 
able  for  operations  at  lower  temperatures  because  of  the 
substantially  lower  ohmic  resistance  of  the  electrolyte.  By 
lowering  the  operating  temperature  of  SOFCs  below  800  °C, 
conventional  metal  alloys  can  be  used  for  interconnectors 
or  auxiliary  components  with  which  a  high  mechanical  re¬ 
liability  of  a  cell-stack  and  lower  manufacturing  costs  can 
be  achieved.  Thus,  recently,  anode-supported  SOFCs  are  re¬ 
ceiving  much  interest  [1-5].  The  anode-supported  cells  are 
generally  fabricated  by  co-firing  the  thin  electrolyte  on  the 
anode  substrate  at  1400-1500 °C  [1-5].  Ni/YSZ  and  YSZ 
are  used  as  the  anode  substrate  and  the  electrolyte,  respec¬ 
tively.  The  thermal  expansion  coefficients  (TEC)  of  Ni/YSZ 
and  YSZ  do  not  match  well,  and  as  a  result  the  mismatch  in¬ 
duces  a  large  residual  stress  in  the  cell  at  room  temperature. 
The  thermal  stress  may  induce  micro-cracks  or  a  delamina¬ 
tion  of  the  electrolyte  at  the  interface  under  a  thermal  cycle. 
These  damages  of  the  cell  decrease  the  cell  performance,  and 
sometimes,  cause  a  destruction  of  the  cell.  Thus  it  is  signif¬ 
icant  to  reduce  the  thermal  stress  under  the  thermal  cycles. 

We  have  been  studying  the  cell  performance  and  the  ther¬ 
mal  stresses  in  the  cell  using  computer  simulation  for  the 
self-supporting  type  planar  SOFCs  [6-8].  Recently  we  have 
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also  carried  out  numerical  analysis  on  cell  performances 
[9]  and  thermal  stresses  on  the  anode-supported  cell.  How¬ 
ever,  for  the  simulation  of  the  thermal  stress,  precise  calcu¬ 
lations  are  difficult  to  perform  because  the  residual  stresses 
in  the  cell  at  room  temperature  are  not  known  well.  Hence 
estimating  the  residual  stresses  in  the  electrolyte  of  the 
anode-supported  cells  at  room  temperature  is  useful  for  the 
calculation  of  the  thermal  stresses  in  the  anode-supported 
cell  under  cell  operations. 

In  this  paper,  we  report  the  evaluation  of  the  residual 
stresses  in  the  electrolyte  of  the  anode-supported  cell  at 
room  temperature.  The  X-ray  diffraction  method  was  used 
to  measure  the  residual  stresses  in  the  electrolyte  of  the 
anode-supported  cell,  and  a  synchrotron  radiation  was  used 
as  an  excellent  X-ray  source,  which  enables  the  estimation  of 
the  residual  stresses  in  the  electrolyte  with  a  high  accuracy. 
We  also  carried  out  a  numerical  calculation  for  the  residual 
stresses  in  the  electrolyte  of  the  anode-supported  cell,  and 
made  clear  that  the  calculated  stresses  were  close  to  those 
experimentally  measured. 


2.  Experimental 

We  measured  the  X-ray  diffraction  pattern  for  the  elec¬ 
trolyte  of  the  anode-supported  cell,  and  estimated  the  resid¬ 
ual  stresses  by  the  sin2  \[r  method  [10-14],  When  a  sample 
is  sustaining  stresses,  the  interplanar  spacing  d  between  the 
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specified  diffracting  planes  of  a  crystallite  in  a  microscopic 
grain  can  be  expressed  as 

1  +  V  2  I  V  | 

d  =  — — — do(J  sin  i/r  +  {1  -  —  (cri  +er2)jr/o  (1) 

where  do  is  the  interplanar  spacing  under  a  stress-free  con¬ 
dition,  ij/  the  tilting  angle  of  the  specified  diffracting  plane 
in  the  crystallite  from  the  normal  to  the  sample  surface,  v 
the  Poisson’s  ratio,  E  the  Young’s  modulus,  a  the  stress  in 
the  grain,  and  cs\  and  01  are  the  principal  stresses  [10-13], 
Modifying  Eq.  (1),  a  is  expressed  by 

E  1  dd 

ct  = - (2) 

( 1  +  v)  do  9sin“  i fr 

To  derive  Eq.  (2),  it  is  assumed  that  the  sample  is  isotropic 
and  elastic.  In  addition,  there  is  no  stress  perpendicular  to 
the  sample  surface,  because  only  the  stresses  near  the  sample 
surface  are  measured  using  the  X-ray  diffraction  method. 

In  this  stress  measurement  we  employed  the  iso-inclination 
method  and  the  fixed  i/'p  method.  The  configuration  for 
the  diffraction  measurements  and  the  stress  analyses  is  in¬ 
dicated  in  Fig.  1.  In  Fig.  1,  t/t0  is  the  angle  between  the 
normal  to  the  sample  surface  and  the  incident  beam  and  26 
is  the  scattering  angle  of  the  specified  diffracting  plane  in 
the  crystallite.  In  the  fixed  i/'p  method,  26  scan  is  carried 
out  at  a  fixed  i/fo-  We  determine  +ij/  and  — i//  as  follows 
(+i/r;  the  sample  normal  is  tilting  to  the  diffraction  beam, 
— i j/\  the  sample  normal  is  tilting  to  the  incident  beam).  To 
eliminate  the  effect  of  an  orientation  of  the  crystalline  or 
an  irregularity  of  the  grain  on  the  stress  measurements,  the 
sample  was  fluctuated  by  ±1°  around  the  fixed  t//o  with  an 
enough  high  speed  compared  with  the  scanning  speed  of  26. 

A  synchrotron  radiation  was  used  as  an  X-ray  source.  The 
synchrotron  radiation  is  an  ultra-bright,  highly  directional, 
and  collimated  light.  In  addition,  a  wide  spectrum  of  wave¬ 
lengths,  stretching  from  infrared  to  X-rays  can  be  used.  Us¬ 
ing  the  synchrotron  radiation  as  the  beam  source  in  the  resid¬ 
ual  stress  measurements  has  the  following  advantages  [14]: 


1.  The  beam  with  an  arbitrary  wavelength  is  available  and 
thus  the  stress  measurement  for  an  arbitrary  diffracting 
plane  with  a  high  accuracy  is  possible. 

2.  The  beam  is  the  monochromatic  light  and  thus  it  is  not 
necessary  to  consider  the  effect  of  Ka iiC(2  doublet. 

The  experiment  was  performed  at  BL09XU  line  at 
SPring-8,  which  is  a  standard  X-ray  undulator  beam  line. 
For  the  stress  measurements,  we  selected  an  X-ray  energy 
of  8.05  keV  (X  —  0.1541  nm),  which  corresponds  to  the  Cu 
Ka  line,  to  compare  the  present  results  with  those  mea¬ 
sured  using  a  commercial  X-ray  diffraction  apparatus.  The 
attenuation  length  of  the  X-ray  of  X  =  0.1541  nm  in  ZrCU 
at  a  fixed  angle  of  90°  is  calculated  to  be  around  15  p,m. 
The  typical  thickness  of  the  electrolyte  was  around  30  p,m, 
which  is  almost  triple  the  attenuation  length  of  the  X-ray 
in  ZrCF,  and  thus,  the  contribution  of  the  anode  to  the 
diffraction  beam  was  negligible. 

The  beam  size  at  the  sample  point  was  adjusted  to 
1  mm  ( H)  x  1  mm  ( V)  using  W  slits.  Regarding  the  diffract¬ 
ing  plane  used  for  the  stress  measurements,  we  selected 
(5  3  1)  plane  of  8YSZ. 

The  anode-supported  cells  used  for  the  stress  measure¬ 
ments  were  fabricated  by  co-firing  both  the  electrolyte 
and  the  anode  at  1500 °C  [15].  The  8YSZ  electrolyte  was 
first  screen  printed,  or  dip-coated  on  the  green  NiO/YSZ 
substrate.  To  avoid  generations  of  cracks  and  micro-pores 
in  the  electrolyte,  the  particle  sizes  of  the  raw  powders 
for  NiO/YSZ  substrate  and  YSZ  electrolyte  were  selected 
carefully.  After  co-firing,  the  anode-supported  cell  was 
cooled  down  to  room  temperature.  The  typical  cell  size 
was  40  mm  x  40  mm  x  2  mm,  and  the  thickness  of  the 
electrolyte  was  about  20-40  |xm  (the  image  is  shown  in 
Fig.  2).  As  compared  to  the  anode  substrate,  the  electrolyte 
is  thin  enough  to  assume  that  the  direction  of  the  main 
stress  in  the  electrolyte  is  parallel  to  the  film  surface.  To 
investigate  the  stress  distribution  in  the  anode-supported 
cell,  two  10  mm  x  10  mm  pieces  were  cut  from  the  center 
and  the  corner  parts  of  the  sample.  In  this  measurement,  we 


Surface  normal 


Fig.  1.  Schematic  diagram  of  the  X-ray  stress  measurement;  iso-inclination  method  and  fixed  i/tq  method. 
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Electrolyte  (8YSZ) 


Fig.  2.  Schematic  view  of  the  anode- supported  cell. 


Table  1 


List  of  the  properties  of  the  samples  used  for  the  stress  measurement 


ID 

Electrolyte 

material 

Manufacturing 

process 

Post 

treatment 

Anode 

state 

Position 

1 

8YSZ 

Screen-print 

As-sintered 

Oxidized 

Center 

2 

8YSZ 

Screen-print 

As-sintered 

Oxidized 

Comer 

3 

8YSZ 

Screen-print 

Flattened 

Oxidized 

Center 

4 

8YSZ 

Screen-print 

Flattened 

Oxidized 

Corner 

5 

8YSZ 

Screen-print 

Flattened 

Reduced 

Center 

6 

8YSZ 

Screen-print 

Flattened 

Reduced 

Comer 

7 

8YSZ 

Dip-coat 

Flattened 

Oxidized 

Comer 

measured  the  residual  stresses  for  seven  different  samples. 
The  characteristics  of  the  samples  are  listed  in  Table  1. 
Samples  1  and  2  were  cut  from  an  as-grown  cell.  Usually, 
the  as-grown  cell  warped  toward  the  electrolyte  side  as 
shown  in  Fig.  3.  If  the  cells  warp,  it  is  difficult  to  contact 
each  cell  well,  and  the  electric  connection  between  cells  be¬ 
comes  poor.  Hence,  we  reformed  the  warp  and  flattened  the 
cell  using  a  distributed  load  on  the  electrolyte  at  1400  °C. 


Ni 


- ^  y*Lx 

3.524  A 

Fig.  4.  Comparison  of  the  lattice  structure  and  the  relevant  lattice  constant 
between  NiO  and  Ni. 

Samples  3  and  4  were  cut  from  the  flattened  cell.  Before 
an  operation  of  the  cell,  the  anode  part  of  the  cell  has  to  be 
exposed  to  a  reducing  atmosphere  to  reduce  NiO  to  Ni  in 
the  anode.  Samples  5  and  6  were  the  reduced  samples.  The 
lattice  constant  of  Ni  is  smaller  by  20%  than  that  of  NiO  as 
illustrated  in  Fig.  4.  As  a  result,  the  reduction  of  the  anode 


NiO 


4.168  A 


Fig.  3.  Side  view  of  the  warped  anode-supported  cell.  The  electrolyte  is  located  at  the  upper  side. 
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Table  2 


List  of  the  mechanical  properties  and  cell  size  used  for  the  stress  simulation 


Young’s 

modulus 

(GPa) 

Poisson’s 

ratio 

TEC  (K_1) 

Thickness  (mm) 

Electrolyte 

206“ 

0.3b 

10.56  x  10~6c 

40  x  40  x  0.03 

Anode 

96d 

0.3b 

12.22  x  10~6c 

40  x  40  x  2 

“Taken  from  Ref.  [16]. 

b  Assumed  value. 

“Assumed  from  the  temperature  dependence  of  TEC  below  1000 °C 

[8], 

d  From  4-point  bending  tests  at  room  temperature. 

causes  a  shrinkage  of  the  anode,  and  may  change  the  resid¬ 
ual  stress  in  the  electrolyte.  For  Samples  1-6,  the  electrolyte 
is  screen  printed  on  the  surface  of  the  anode  substrate  first 
and  then  co-fired.  Only  for  Sample  7,  the  electrolyte  was 
dip-coated  on  the  surface  of  the  anode  substrate. 

Numerical  calculations  for  the  thermal  stresses  in  the 
anode-supported  cells  were  carried  out  using  the  finite  ele¬ 
ment  program  “ABAQUS”  (Hibbit,  Karlsson  and  Sorensen 
Inc.).  We  modeled  the  geometry  of  the  cell  of  the  elec¬ 
trolyte/anode  bi-layer,  and  calculated  the  residual  thermal 
stresses  at  room  temperature.  In  the  calculation,  we  assumed 
that  both  the  electrolyte  and  the  anode  were  constrained  each 
other  below  1400  °C,  and  the  origin  of  the  residual  stresses 
in  the  cell  was  only  the  mismatch  of  TEC  between  the  elec¬ 
trolyte  and  the  anode.  The  mechanical  properties  and  the 
cell  size  used  for  the  stress  calculation  are  listed  in  Table  2. 

3.  Results  and  discussion 

Fig.  5  shows  the  result  of  0-26  scan  for  the  electrolyte  of 
the  Sample  3.  First  the  sample  was  set  so  that  the  sample 
surface  was  parallel  to  the  incident  beam,  and  then  the  inci¬ 


dent  beam  was  spotted  at  the  center  part  of  the  sample.  All 
the  peaks  in  the  diffraction  pattern  can  be  assigned  to  the 
cubic  YSZ  structure.  There  are  no  diffraction  peaks  from 
the  anode  substrate. 

Fig.  6  displays  the  change  of  the  diffraction  peak  of  YSZ 
(5  3  1)  plane  with  the  change  of  i]j  measured  for  the  Sam¬ 
ple  3.  With  increasing  i/f,  the  peak  position  shifts  to  the 
higher  angle,  indicating  that  the  inter-planar  spacing  cl  de¬ 
creases  and  the  stress  type  is  compressive.  To  determine 
the  center  of  the  peak,  the  peak  profiles  are  fitted  with  the 
Gaussian-function,  and,  thus  produced  af-sin2  xji  diagram  is 
shown  in  Fig.  7.  The  open  circles  and  open  squares  are  for  the 
positive  and  the  negative  i/r,  respectively,  and  the  solid  line 
is  the  linear  least-squares  fit.  Although  the  data  are  slightly 
scattered,  d  decreases  linearly  with  the  increase  of  sin2  i jr. 
From  the  slope  of  the  c/-sin2  i jr  diagram,  the  residual  stress 
can  be  estimated  using  Eq.  (2).  In  the  X-ray  stress  mea¬ 
surements,  the  stress  in  a  sample  can  be  estimated  using  the 
diffraction  from  the  specified  diffracting  plane  of  crystallites 
which  orients  to  the  specified  direction.  Therefore  the  elastic 
constant  E  and  v  in  Eq.  (2)  is  different  from  those  measured 
by  a  mechanical  method,  and,  this  elastic  constant  should  be 
measured  using  the  X-ray  diffraction  method.  In  this  study, 
however,  it  was  difficult  to  measure  E  and  v  for  the  prac¬ 
tical  cell  by  the  X-ray  diffraction  method.  Since  there  are 
no  data  for  the  X-ray  elastic  constant  measured  using  YSZ 
(5  3  1)  plane,  we  calculated  the  elastic  constant  assuming  the 
Voight  model  [12,17]  using  the  reported  elastic  constant  for 
cubic  YSZ  [18].  The  reported  elastic  constant  for  cubic  YSZ 
and  the  calculated  £7(1  +  v)  are  listed  in  Table  3.  For  the 
Sample  3,  the  estimated  residual  stress  in  the  electrolyte  was 
a  compressive  stress  of  633  MPa.  For  all  the  other  samples, 
the  estimated  residual  stresses  in  the  electrolyte  are  listed  in 
Table  4.  From  the  present  results,  the  following  information 
was  obtained. 


Fig.  5.  X-ray  diffraction  pattern  for  the  electrolyte  of  the  anode- supported  cell. 
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20  (deg) 


Fig.  6.  Change  of  the  diffraction  peak  of  YSZ  (5  3  0)  plane  with  the  change  of  t/r  for  the  Sample  3. 


Fig.  7.  rf-sin2  <!/  diagram  measured  for  the  Sample  3;  open  circles  and  open  squares  are  for  the  positive  and  the  negative  i/r,  respectively,  solid  line  is 
the  linear  least-squares  fit.  d  was  determined  from  the  diffraction  peak  of  YSZ  (5  3  0)  plane. 


(a)  By  comparing  1  and  2,  the  residual  stress  was  larger  at 
the  center  part  of  the  cell.  The  origin  of  the  difference 
was  the  macro  deformation  of  the  cell  (i.e.  the  warp  of 
the  cell). 

(b)  By  the  flattening  treatment,  the  residual  stress  at  the 
center  part  of  the  cell  was  lowered,  and  the  distribution 
of  the  residual  stress  in  the  electrolyte  became  small. 

Table  3 


List  of  the  stiffness  of  the 

elastic  constant 

cubic  YSZ  used  for  the 

calculation  of  the 

Cn  (GPa) 

C12  (GPa) 

C44  (GPa) 

E/(  1  +  v) 

(x  10”4  GPa-1) 

Electrolyte 

400a 

110“ 

55a 

182 

This  means  that  the  electrolyte  is  in  a  plastic  state  at 
1400 °C. 

(c)  The  effect  of  the  reducing  treatment  of  the  anode  on  the 
residual  stress  is  small.  This  may  come  from  that  the 
Young’s  modulus  of  the  anode  became  small  owing  to 
the  reduction.  That  is,  the  effects  of  the  softening  of  the 
anode  may  just  cancel  the  effect  of  the  expansion  on  the 
residual  stress. 


Table  4 


List  of  the  results  of  the  simulation  for  the  residual  stress  in 

the  electrolyte 

Sample  ID 

1  2 

3 

4 

5 

6 

7 

Calculated  stress 
(MPa) 

774  672 

633 

677 

658 

678 

676 

Taken  from  Ref.  [18]. 
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Fig.  8.  Simulated  stress  distribution  in  the  electrolyte  at  room  temperature. 


(d)  The  residual  stress  is  independent  of  the  fabrication 
method  of  the  electrolyte.  This  is  reasonable  because  the 
origin  of  the  residual  stress  may  be  only  the  mismatch 
of  TEC  between  the  electrolyte  and  the  anode. 

Numerical  simulations  for  the  residual  stresses  in  the 
cell  at  room  temperature  were  also  carried  out  to  check 
the  validity  of  the  stress  measurement.  The  sample  geom¬ 
etry  used  for  the  simulation  was  the  same  as  the  practical 
cell  used  for  the  X-ray  stress  measurements.  Fig.  8  shows 
the  simulated  stress  distribution  in  the  electrolyte  for  the 
anode-supported  cell.  This  simulation  corresponds  to  the 
flattened  sample.  Except  for  the  edge  part,  the  stresses 
are  almost  homogeneous  values  of  640-700  MPa  over  the 
electrolyte  plane.  Near  the  corner  part  the  stress  is  slightly 
larger  than  that  at  the  center  part,  which  is  qualitatively 
consistent  with  the  result  of  the  X-ray  stress  measure¬ 
ments.  In  this  calculation,  the  electrolyte  and  the  anode 
is  assumed  to  be  constrained  to  each  other  at  1400  °C. 
However,  the  temperature  at  which  the  electrolyte  and  the 
anode  are  constrained  is  not  clear,  and  the  assumption 
of  the  temperature  affects  the  results  of  the  simulation 
as  shown  in  Fig.  9.  With  the  change  of  the  temperature 
at  which  the  electrolyte  and  the  anode  are  constrained, 
the  calculated  stress  at  the  center  part  varies  from  550  to 
670  MPa.  From  the  following  fact,  the  temperature  at  which 
the  electrolyte  and  the  anode  are  constrained  seems  to  be 
between  1250  and  1400  °C.  The  cell  could  be  flattened  at 
1400  °C,  and,  the  present  stress  measurement  revealed  that 
the  flattening  treatment  did  not  affect  the  residual  stress, 
suggesting  that  the  electrolyte  and  the  anode  may  be  con¬ 
strained  below  1400  °C.  In  contrast,  it  was  reported  that 
the  sintering  of  YSZ  does  not  proceed  well  below  1250  °C 
[19].  When  the  temperature  at  which  the  electrolyte  and 
the  anode  are  constrained  is  assumed  to  be  1400  °C,  the 
calculated  residual  stress  is  close  to  the  measured  value. 
Thus,  in  the  following  stress  calculations  we  assumed 


Fig.  9.  Effect  of  the  temperature  at  which  the  electrolyte  and  the  anode 
are  constrained  on  the  calculated  residual  stress  in  the  electrolyte. 


that  the  electrolyte  and  the  anode  would  be  constrained  at 
1400 °C. 

Here  we  will  focus  on  the  effect  of  the  cut  of  the  sample 
on  the  residual  stress.  In  the  stress  measurements,  the  cell 
was  cut  into  the  10  mm  x  10  mm  sample.  This  came  from 
the  restriction  of  the  size  of  the  sample  stage.  It  should  be 
taken  into  account  that  the  cut  of  the  cell  would  change  the 
residual  stress  in  the  sample.  Hence  using  the  stress  simu¬ 
lation  we  studied  the  effect  of  the  sample  cut  on  the  resid¬ 
ual  stress.  In  the  present  simulation,  the  thicknesses  of  the 
electrolyte  and  the  anode  were  fixed  to  30  |xm  and  2  mm, 
respectively,  and  the  residual  stress  was  calculated  with  var¬ 
ious  plane  sizes.  The  result  of  the  simulation  is  displayed 
in  Fig.  10.  When  the  cell  size  is  more  than  5  mm,  the  cal¬ 
culated  stress  is  almost  independent  of  the  cell  size.  In  the 
practical  stress  measurements,  the  sizes  of  the  samples  were 
larger  than  5  mm,  and,  if  macro  deformations  of  the  sample 


Cell  size  (mm) 


Fig.  10.  Effect  of  the  plane  sample  size  on  the  calculated  residual  stress 
in  the  electrolyte.  The  thicknesses  of  the  electrolyte  and  the  anode  are 
fixed  to  30  |xm  and  2  mm,  respectively. 
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Fig.  1 1 .  Effect  of  the  thicknesses  of  the  electrolyte  and  the  anode  on  the  calculated  stress  in  the  electrolyte.  The  plane  sample  size  is  fixed  to  40  mm  x  40  mm. 


(e.g.  the  warp  of  the  cell)  hardly  affect  the  residual  stress, 
the  effect  of  the  sample  cut  on  the  residual  stress  may  be 
negligible. 

Next,  we  will  comment  on  the  effect  of  the  thickness  of 
the  electrolyte  and  the  anode  on  the  residual  stresses.  The 
residual  stresses  in  the  electrolyte  were  calculated  for  the 
difference  combination  of  the  thicknesses  of  the  anode  and 
the  electrolyte,  and  the  results  were  displayed  in  Fig.  11. 
The  combination  of  the  thin  anode  and  the  thick  electrolyte 
simply  decreases  the  residual  stress  in  the  electrolyte.  On  the 
other  hand,  for  a  practical  cell,  the  thick  electrolyte  causes 
an  increase  of  the  electric  resistance  in  the  electrolyte.  Thus 
the  optimal  thickness  of  the  electrolyte  has  to  be  selected 
with  consideration  of  both  the  electric  resistance  and  the 
residual  stress. 

From  the  X-ray  stress  measurements  and  the  stress  sim¬ 
ulation  for  the  electrolyte  of  the  anode-supported  cell,  it 
has  been  made  clear  that  the  residual  stresses  in  the  elec¬ 
trolyte  at  room  temperature  are  compressive  stresses  of 
around  650  MPa.  The  strength  of  YSZ  against  the  com¬ 
pressive  stress  is  more  than  1  GPa  [20]  and  thus  these 
residual  stresses  are  not  so  large  as  to  induce  a  failure  of  the 
electrolyte. 


4.  Conclusion 

The  residual  stresses  in  the  electrolyte  of  the  anode- 
supported  planar  SOFCs  were  estimated  by  experimental 
and  numerical  methods.  The  stress  measurements  were 
performed  by  the  X-ray  stress  measurement,  and  the  nu¬ 
merical  calculations  of  the  stress  were  carried  out  by  the 
finite  element  method.  The  estimated  residual  stresses 
in  the  electrolyte  were  compressive  stresses  of  around 
650  MPa.  The  effects  of  the  cell  geometry,  the  flattening 


treatment,  the  fabrication  method  of  the  electrolyte,  and 
the  reduction  of  the  anode  on  the  residual  stress  were  also 
studied. 
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